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Введение. Статья посвящена теоретическому изучению 
переноса ионов бинарной соли с учетом реакции диссоци- 
ации (рекомбинации) воды. Цели исследования: построе- 
ние математической модели; разработка алгоритма чис- 
ленного решения краевой задачи, соответствующей мате- 
матической модели; развитие теории подобия, включая 
переход к безразмерному виду с использованием харак- 
терных величин; определение физического смысла триви- 
альных критериев подобия; нахождение нетривиальных 
критериев подобия; построение и анализ вольтамперной 
характеристики (ВАХ). 

Материалы и методы. При теоретическом исследовании 
и численном анализе переноса ионов бинарной соли учи- 
тывается реакция диссоциации (рекомбинации) воды. При 
этом используются уравнение теплопроводности и мате- 
матическая модель электродиффузии одновременно четы- 
рех видов ионов (двух ионов соли, а также ионов Н* и 
ОН”) в диффузионном слое электромембранных систем с 
идеально селективной мембраной. 

Для дифференциальных уравнений первого порядка ста- 
вится сингулярно возмущенная краевая задача. В уравне- 
нии для напряженности электрического поля правая часть 
не зависит от напряженности. При численном решении 
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ваются недостающие краевые условия для напряженности 
электрического поля. 

Результаты исследования. Разработаны новая математи- 
ческая модель, алгоритм численного решения краевой 
задачи, программное обеспечение. Проведен численный 
анализ и определены фундаментальные закономерности 
переноса ионов соли с учетом реакции диссоциации (ре- 
комбинации) молекул воды, температурных эффектов и 
Джоулева разогрева. Построена и проанализирована ВАХ. 
Обсуждение и заключение. Рассмотрен перенос ионов 
бинарной соли через диффузионный слой у катионооб- 
менной мембраны. Предложена математическая модель 
названного процесса, которая учитывает температурные 
эффекты, обусловленные реакциями диссоциации (реком- 
бинации) молекул воды и Джоулевым нагревом в раство- 
ре. Установлены основные закономерности переноса 
ионов соли с учетом реакции диссоциации (рекомбина- 
ции) молекул воды и температурных эффектов. Темпера- 
турные эффекты от реакции диссоциации (рекомбинации) 
и Джоулева разогрева в ОЭН практически незаметны (ис- 
ключение — область рекомбинации, ОР). Джоулев нагрев 
в области пространственного заряда (ОПЗ) на два порядка 
больше охлаждающего эффекта реакции диссоциации 
воды. При рекомбинации в ОР выделяется примерно 
столько же тепла, сколько при Джоулевом нагреве в рас- 
ширенной ОПЗ. Однако из-за малых размеров ОР влияние 
этого тепла незаметно. Значит, можно считать, что есть 
только один источник тепла на межфазной границе в ОПЗ, 
который благодаря своему заметному размеру обусловли- 
вает значительное повышение температуры во всем диф- 
фузионном слое. Отсюда следует, что возможно возник- 
новение и развитие гравитационной конвекции. Общие 
выводы, следующие из полученных результатов, открыва- 
ет возможность интенсификации процесса переноса ионов 
соли в электродиализных аппаратах. 


Ключевые слова: мембранные системы, ионообменная 
мембрана, пространственный заряд, расширенная область 
пространственного заряда, диссоциация (рекомбинация) 
молекул воды. 
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Введение. Интенсивный токовый режим предполагает использование токов, в несколько раз превыша- 


ющих предельный электродиффузионный. В этом случае в мембранных системах возникают вторичные (или 


сопряженные) явления концентрационной поляризации: 


— пространственный электрический заряд занимает макроскопическую область, сопоставимую с тол- 


щиной диффузионного слоя; 


— меняется рН раствора [1-3], что объясняется расщеплением воды; 


— формируются микроконвективные течения и пр. 


Учет влияния реакции диссоциации (рекомбинации) молекул воды важен для понимания процессов пе- 


реноса в электромембранных системах, поскольку ряд авторов считает, что появление новых носителей заряда 


Н* и ОН` может обусловить уменьшение или даже исчезновение пространственного заряда, являющегося ос- 


новой других механизмов переноса (например, электроконвекции). Реакция диссоциации молекул воды — эн- 
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дотермическая, а реакция рекомбинации — экзотермическая. Центры реакционных областей разделены в про- 
странстве. Это вызывает неравномерное распределение температуры и, следовательно, влияет на все физиче- 
ские свойства воды: коэффициент диссоциации, теплоемкость и пр. Кроме того, неравномерное распределение 
температуры может вызвать гравитационную конвекцию. Таким образом, исследование температурных эффек- 
тов, связанных с реакциями диссоциации и рекомбинации молекул воды, является актуальной задачей [4]. В 
данной работе ограничимся учетом влияния неравномерного распределения температуры на коэффициент дис- 
социации молекул воды. 

Материалы и методы 

1. Математическая модель 

1.1. Система уравнений 

Рассмотрим одномерный стационарный процесс переноса для однозарядных ионов соли в диффузион- 
ном слое Нернста. Учтем при этом реакции диссоциации (рекомбинации) воды и температурные эффекты, свя- 
занные с Джоулевым разогревом раствора. Такой случай описывается системой уравнений [5]: 
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к +Е-1 + В(СзС4 — Ки(Т)) = 0. (7) 
Здесь В. = В» = 0,В; = В. = КаСн»о — К„СзСа = К, (К, — СзС4), Ки(Т) = К ‚еР(Т-То), + = 1,2 — индексы ионов 
соли; { = Зи [= 4 — соответственно индексы ионов водорода Н* и гидроксила ОН”; Е = - — напряжен- 
ность электрического поля; С;, д, О; — соответственно концентрация, поток, коэффициент диффузии 1-го иона; 


Ка(Г) — константа скорости диссоциации воды; К.. — константа скорости рекомбинации ионов водорода Н* и 

гидроксила ОН”; К„— константа равновесия; =. — диэлектрическая проницаемость раствора; Е — число Фа- 

радея; Г — общая плотность тока; К — коэффициент теплопроводности; ро — плотность раствора; с, — удель- 

ная теплоемкость раствора; @ = Ё [Г — плотность источников тепла, связанного с прохождением электриче- 
Касн2о 

ского тока через раствор (Джоулев разогрев раствора); В = акК.., О = В (с; 4 — ие) = “К, 16: 5. (Г)) — 


плотность источников и стоков тепла, связанного с реакциями рекомбинации и диссоциации молекул воды. 
Общая плотность тока Г равна: [ = [$ + [,,. 
5 =Е (1 -]2), № =Е: Цз — 44). (8) 
Здесь [; — плотность парциального тока по ионам соли; [, — плотность парциального тока по ионам 
ВОДЫ. 
Из (6)-(8) следует: и = соп5$®, }› = соп5®, [ =ЁЕ` (1 -]2) = соп5$, [, = с0п5$6, но }з = с0П$Е, 4 = 
К 
СОП5Е, а? = ти коэффициент температуропроводности раствора. 
обр 


1.2. Краевые условия 

Будем считать, что х = 0 соответствует глубине раствора, где выполняется условие электронейтраль- 
ности С: — С. + С. — С, = 0, ах = 6 соответствует условной межфазной границе «раствор — мембрана». Та- 
ким образом, 5 — толщина диффузионного слоя Нернста. 

В этой работе для определенности рассматривается катионообменная мембрана (КОМ), которая для 


<“ „(14 
простоты будет предполагаться идеально селективной . = 0). 





Рис. 1. Схема диффузионного слоя у катионообменной мембраны 
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1) Краевые условия в глубине раствора (х = 0) 
Положим, С1(0) = Сло, С2(0) = Соо, Сз(0) = Сзо, С4(0) = —— "> причем Сло + Сзо = С2о, Т(6) = То. 


2) Краевые условия на межфазной границе И (х=б) 

Полагаем, С, (6) = С\т, С3(6) = Сзт. При этом сумма Ст + Сът сопоставима с табличным значением 
обменной емкости катионообменной мембраны. Плотность потока (5) = мт ионов ОН `` определяется ката- 
литической реакцией диссоциации молекул воды в тонком поверхностном слое катионообменной мембраны. 
Для температуры используется условие баланса потоков в точке х = 6 [1]: 

НКТ = 6(5)6 + 0(6)5 +49 

2. Переход к безразмерному ви 

2.1. Характерные величины и формулы перехода к безразмерному виду 

В качестве характерной концентрации Су возьмем концентрацию катионов соли С10о. Реальные значения 
С1о находятся в диапазоне 10-?М -— 1М. В качестве характерной температуры возьмем комнатную температуру 
То = 298К. В качестве характерной длины рассмотрим толщину диффузионного слоя 6, который имеет порядок 








10 %м. Обозначим Ду = о — коэффициент диффузии электролита (МаС), Фо = = — тепловой потенци- 
1 2. 
ал [5]. 
о к безразмерным величинам: 
хо = И фе № Е = —6 ев, 19) = — 2% = > 1% = —— Е =”,  ^ — _- _ — о = 
0 Фо ах пр пр” тр Сто Р:С1о 
ЛЕ — (и) — — бб - 
‚= 1,2,3А, т =- 1, = 
Е т 1 пр 5 
(м) —_ КиСГ) _ (то — —_ Км —_ 
Тогда Ки” (Г) = и = уе (т 1). где у@2 2 р = ЬТь. 
О АИ ты 
Положим, & С оЕ252° т 


2.2. Переход к безразмерному виду в уравнениях 
Ориентируясь на [6], перейдем к безразмерному виду в системе уравнений (1)—(6): 


(и) 
аб 1 _ (м) .(и) 
= С: Е д“, (9) 
(и) 
а(> ^ _ _ (и) (и) _ и 
о (10) 
0 
асз— (и) (и) _ (4) 
и Е (11) 
(и) 
ао _ _ „(м . (и) 
пе о. (13) 
0 р@д ры ыы (Рева) —_ 6. (14) 
а _ (м) _ „ад, 0) _ пы) 
В о (15) 


Запишем общую плотность тока в виде 1%) + 1 = 1%. 


Здесь 


16) № вы в” 


(1) 


= 


= 0 и, следовательно, 
1%) — 2“ (и) 


причем для идеально селективной мембраны ]> 








Л 
При переходе к безразмерному виду уравнение теплопроводности (7) приобретает вид: 
р (и) (ти) — 
ре ЕО 199 + ВО (СС — удев" 1-9) = 0. (16) 
Здесь К(Э = 60 Куб” С1о 
РоС1оВ РовТо ° 


2.3. Физический смысл и оценка величин тривиальных критериев подобия 

Параметры =), В р а, ко, К 2, В05, Ь@®) являются тривиальными критериями подобия. Опре- 
делим их физический смысл и значения при характерных значениях входных параметров задачи. 

1. 5% = 15. 

2. Параметр у представляет собой квадрат отношения концентрации ионов Н* или ОН” в нейтраль- 
ном растворе (ФН = 7) к концентрации электролита. Его значения меняются в пределах от 10—10 до 10-18, ион 
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можег считаться малым параметром. 





2 

ВТо= 1 ВТЕ 

3. = = о. =” (4) ‚ где [1 = | > — дебаевская длина. Таким образом, параметр =@9 является 
10 10 


удвоенным квадратом отношения дебаевской длины [1 к ширине диффузионного слоя б. Значения = меняют- 
ся в пределах от 10-11 до 1077, и этот параметр можно считать малым. 





2 

ВТоЕу К ВТо=у" 2КуС [ р 

Е бы Сы (№) ‚ Здесь [.. = `_. Данная величина пропорциональная ширине 
Е? Ро 2б10Е? Ро Е бо 


области рекомбинации (ОР). Таким образом, параметр а@ — это квадрат отношения дебаевской длины к ши- 


рине ОР [., . Параметр а = 14 ив рамках данной задачи является универсальной постоянной, поскольку не 


[ 
зависит от входящих параметров задачи. Таким образом: _ = \/а 9 = 3,74, т.е. ширина ОР составляет при- 
т 


мерно 37 % от ширины квазиравновесной области пространственного заряда (ОПЗ, дебаевская длина). 
2.4. Нетривиальные критерии подобия, физический смысл и оценка величин 
В уравнении (16) безразмерные коэффициенты К@9 и В© являются тривиальными критериями подо- 
бия. Любое соотношение между тривиальными критериями называется нетривиальным критерием подобия. 
Ниже выразим К и В через другие известные критерии подобия. 
К а’роср _ а? Роср 2 


Рос © 
1) Рассмотрим КО = = ОР. ОР =. ОР Здесь — — число Льюиса — Семенова. Введем 
Роб1оК  Боб10о По С1оВ СтоК Бо 





Со 1 
в рассмотрение безразмерное число и = ЕЕ ‚ Тогда К@9 =. -. 
оСр и 


Число и обозначает относительную удельную объемную теплоемкость, изменяется в пределах от 10-8 
до 10“, и его можно рассматривать как малый параметр. 








ак!-52 С а5? 2К..С а 52 52 а 
2) Рассмотрим параметр В = = "^_1® = —5— = Атп--. Здесь Ати = — = 22,8 — 
Ровто — 2ЕТо Бо 2ЕТо о 212 ВТо 
те10 


безразмерная величина, являющаяся тепловым аналогом критерия Аррениуса. Она характеризует чувствитель- 


2 
> 1 
ность количества тепла протекания химическои реакции к изменению температуры. Так как (>) — 


нада 2 ово ди 


Таким образом, уравнение теплопроводности (16) имеет следующий безразмерный вид: 








ат 
&Ёе 
3) Вычисления показывают, что 
= К» В _ КиЕЕ? 6" В _ Ку 62 
Сто РоСбр ВТо=у РоСр РоСрТоЕт 


Ки„Е262 
РоСрТо т. 


терий подобия, который показывает соотношение 


2 
1 и е 
Поскольку #09 = 2 (4) ‚ то М = = 2,5 . 10-3. Следовательно, получаем нетривиальный кри- 


и. 


ун = Па. 
Аналогично 
ЕН = Пал 
Здесь 
ВАТ —14 
ПЕ 30.10 17 
а,1 росрЕ? 52 , ( ) 
— 2 
= Ма 2Е > (18) 
2464 К 
М2 = — 
В2Т2 Е? 


Заметим, что п2 = Па1 Ма. 

2.5. Безразмерный вид системы уравнений 

Таким образом, безразмерная система уравнений, описывающая модель, имеет вид (индекс и опущен 
для простоты записи): 


ас ‚ ас ‚ ас у 

к = СЁ — Лох Я Е — рта, 2 СЁ — Л 
а. в ас ; 

ЕБз--)з = а(уе?Т-® — С.С.), = р 


а. — Е 
Е. а = а(уе?-Ъ — С.С,), 52. = С, — С> + С3 — Са. 


а?т — 
Ге > + ЕНЕТ + цаАтт (СзС. — уе 75) = 0. 
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Система уравнений содержит три малых параметра &, ц и у, связанных соотношениями (17) и (18). Па- 
раметр = входит в уравнения сингулярно, т. е. как коэффициент при производных, а параметр у — регулярно, 
поэтому система уравнений является одновременно сингулярно и регулярно возмущенной. 

2.6. Переход к безразмерному виду в краевых условиях 

Безразмерные краевые условия имеют вид: 


Су|х=о = 1, С>|х=0 = 1 Сзо, Сз|х=о = Сзо, Са|х=о = Со = Г(0) = 1, 
С. х=1 = Слть Сз[=а = Сзт» Ла|х=1 = Ллть 
‚ а Я. 
Ка ГС) = НО4Ба чел. 
Замечание 1. В силу непрерывности справедливо: 
Лу = Эз]зт — Э4Лат = Оз]зо — Вало. 
3. Алгоритм численного решения 
3.1. Переход к краевой задаче для системы уравнений второго порядка 
Из-за наличия малого параметра = в уравнениях (12), (14) и (15) краевая задача (9)—(16) относится к 
сингулярно возмущенным, поэтому для численного решения система неудобна [7]. Кроме того, дифференци- 
альные уравнения имеют первый порядок, а для них ставится краевая задача. В уравнении (15) для напряжен- 
ности электрического поля правая часть не зависит от Е. Этим обусловлены проблемы с устойчивостью при 
численном решении дискретизированной системы уравнений методом Ньютона — Канторовича [8]. Поэтому 


краевую задачу (9)—(16) следует привести к удобному для численного решения виду. Для этого перейдем к си- 
стеме уравнений второго порядка: 








аб: .. @. а2с› _ -а а? _ а — а(уе (Т--с3с.) 
ах? ах (С.Е), ах? ах (С›Е), ах? ах (СзЕ) =Рз ? 
ас, _ а — ат №364) „аЕ _ пр 96 _ 964 
ах? ах (С4Е) =Ра | 2 — Е(С, ы С») тт ах ах 
а?т Ь(Т-1)\\ — 
Ее + НЕТ +: цаАти (С3 С. — уе ) = 0. (19) 


В связи с переходом к системам уравнения второго порядка необходимо дополнить исходные краевые 
условия недостающими краевыми условиями для функции Е. Из краевых условий следует, что в точке х = 0 
выполняется условие электронейтральности С1о — (50 + Сзо — Со = 0. В силу непрерывности с некоторой точ- 
ностью в окрестности точки х = 0 выполняется С; — С. + С. — С, = 0, поэтому Еу может быть вычислено из 
условия электронейтральности. Складывая уравнения для (:, 1 = 1, ...,4 и используя условие электронейтраль- 


+4 (С1(0) — С2(0)). 


_ = (С, — С> + С; — (4) = 1. 


х= 





ности (10 — С50 + Сзо — Со = 0, можно получить соотношение ЁЕ|„-о = 


АЕ 
Из уравнения для Ё прих = 1 имеем: Е 
Хх 


Результаты исследования 
3.2. Результаты численного решения 
Для анализа результатов численного решения наряду с графиками искомых функций С}, Е, }з, ]л, Г рас- 
смотрим графики функций р(х) = С, — С, + С; — С4, р(х) = С.С, — уе Ч-Ъ. Функция р(х) характеризует рас- 
пределение плотности зарядов, а р(х) — отклонение от равновесия реакции диссоциации (рекомбинации). 
Глпе ОтарН:С3*С4-отатта*ехр(Ь*(Т-1)) 


х10 
1.6 
1.5 
1.4 
1.3 
1.2 
1.1 


0.9 
0,8 
0.7 
0.6 
0,5 
0.4 
0,3 
0.2 
0.1 


С3*С4-ватта*ехр(Ь*(Т-1)) 


0 0,1 0,2 0,3 0.4 0,5 0,6 0,7 0,8 0,9 1 
Агтс еп 


а) 
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С1-С2-+С3-С4 


Переп4еп{ уапае Т 


С3*С4-ватта*ехр(Ь*(Т-1)) 
= 





Глпе ОгарБ:С1-С2+С3-С4 


0.006 
0.005 
0.004 
0.003 
0.002 


0.001 


0.001 
—0.002 
0.003 


—0.004 
0 0.1 0.2 0.3 0,4 — 0,5 0.6 0.7 0.8 0,9 
Атс еп 


6) 


Глпе ОтарН:С3*С4-отатта*ехр(Ь*(Т-1)) 


х10- 


—1,2 
—1,4 
—1,6 
—1,8 


0.5 0.81 0,82 0,83 0.84 0,85 0.86 0,87 0,88 0,89 


Агс еп 


в) 
Глпе ОгарВ: Оерепдет{ уапае Т 


1.0028 
1.0026 
1.0024 
1.0022 

1.002 
1.0018 
1.0016 
1.0014 
1.0012 

1.001 
1.0008 
1,0006 
1.0004 
1.0002 





0 0.1 0.2 0.3 0.4 0,5 0.6 0,7 0,8 0,9 1 
Атс еп 


г) 


Рис. 2. График функции р(х) = С, — С + С; - С, (10: 0,93] (а); график функции р(х) = СзС. — уе Т-1(6); 
график функции р(х) = СС. — уе -Ъ в ОПЗ (в); график температуры (2) 


Из рис. 2 а следует, что диффузионный слой делится на три области. Первая — область электроней- 


тральности (ОЭН), где р(х) = 0. Вторая — область пространственного заряда (р(х) > 0). Третья — промежу- 
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точный слой между ними, где плотность распределения зарядов имеет ярко выраженный локальный максимум. 
Сама ОПЗ делится на две части: расширенную и квазиравновесную (она очень мала и примыкает к межфазной 
границе). В квазиравновесной части ОПЗ распределение плотности зарядов подчиняется экспоненциальному 
закону (на рис. 2 не приведено). 

Общий вид графика функции р(х) (рис. 2 6) показывает, что отклонение от равновесия происходит: 

— в области рекомбинации, которая находится внутри ОЭН; 
— в области пространственного заряда (рис. 2 в). 

В первом случае р(х) > 0, то есть рекомбинация преобладает над диссоциацией С3С. » уе? ("-2. Во 
втором случае, наоборот, р(х) < 0, и диссоциация преобладает над рекомбинацией в ОПЗ, а также в промежу- 
точном слое. Температура (рис. 2 г) увеличивается линейно на 0,0028 в безразмерном виде (или 0,0028 : 298 = 
0,83 градуса). Увеличение зависит от общей плотности тока или скачка потенциала в диффузионном слое и 


может составить до десятков градусов. Очевидно, что все члены уравнения (19) по сравнению с первым малы, 
2 


а“т В: 
то есть вторая производная ав = 0. Из этого следует, что температура является линеинои и определяется крае- 
Хх 


выми условиями. 
Глпе ОгарВ: Оереп4деп{ уамае С1 Глпе ОтарВ: Вереп4еп{ уапа е С2 


Глпе ОгарВ: Оереп4деп{ уамае СЗ Глпе ОтарВ: Вереп4еп: уапа ]е С4 


—— С1 (Ма+) 
СС) 
—— С3(Н+) 
—— 64 (ОН-) 





0 0.1 0,2 0.3 0.4 0.5 0.6 0,7 0,8 0.9 1 


Глпе ОгарВ: Оереп4еп{ уамае С1 Глпе ОтарВ: Береп4еп{ уапа е С2 
Глпе ОгарВ: Оереп4еп уамае СЗ Глпе ОтарВ: Ререпдеп{ уапа ]е С4 | 


1,035 

1.03 С1 (Ма+) 
С2 (С1-) 
СЗ (Н+) 

—— 64 (ОН-) 


1.025 
1,02 
1.015 
1,01 
1.005 


0.995 
0,99 
0.985 
0,98 
0.975 
0,97 
0.965 
0.96 
0,955 
0.95 





—0,005 0 0.005 0,01 0,015 0.02 0,025 


0) 
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Рис. 3. Общий вид графиков функций Ст, (>, С5, С (а); графики функций С\ С в окрестности нуля (6); графики функций С3, 
С. в окрестности нуля (в); график функций Ст, С, С3, С. в области пространственного заряда (г) 


Как следует из рис. 3 а, концентрация натрия линейно падает в ОЭН и становится малой в ОПЗ, но она 
выше, чем концентрация остальных ионов (рис. 3 г). Центром ОР является точка пересечения графиков концен- 
трации Ма" и СГ и соответственно ОНТ и Н*. Правее центра концентрация С[` убывает быстрее, чем концен- 
трация Ма". Возникающий из-за этого дефицит отрицательно заряженных ионов компенсируется быстрым ро- 
стом ионов гидроксила, и сохраняется условие электронейтральности (рис. 3 а). 
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Рис. 4. Общий вид графика функции Ё за исключением межфазной границы (а); график функции Ё вблизи ОР (6); 
график функции Ё вблизи межфазной границы (в), вольтамперная характеристика (г) 
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Из данных рис. 4 следует, что в ОЭН напряженность электрического поля быстро увеличивается в де- 
сятки раз (рис. 4 а). Особенно высокая скорость роста напряженности электрического поля отмечается в ОПЗ 


1 2 < 
при достижении величин порядка --. В ОР имеется внутренний погранслой по напряженности электрического 
Е 


поля Ё (рис. 4 6). Это связано с тем, что меняются носители заряда: до ОР заряд определяется ионами Ма*, СГ, 
Н*, а после — ионами Ма", СГ, ОН". Появление «ступеньки» связано с тем, что коэффициент диффузии у 
водорода примерно в два раза больше коэффициента диффузии гидроксила. ВАХ неограниченно растет с уве- 
личением скачка потенциала в запредельном режиме [; > 2 (рис. 4 г). Этот рост происходит значительно мед- 
леннее, чем экспериментально наблюдаемый — значит, наряду с реакцией диссоциации (рекомбинации) воды 
существует и другой механизм сверхпредельного переноса. Общепринято, что таким механизмом является 
электроконвекция [9-11]. В связи с этим возникает проблема оценки и взаимодействия реакции диссоциации 
(рекомбинации) и электроконвекции. 
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Рис. 5. Графики потоков ионов водорода ]з (а) и гидроксила ]^ (6) 

Из данных рис. 5 видно, что поток ионов водорода положительный. Он постоянен слева от ОР и прак- 
тически равен нулю справа от ОР вплоть до ОПЗ. В ОПЗ поток ионов водорода линейно увеличивается в связи 
с интенсивной диссоциацией воды. Поток ионов гидроксила левее ОР практически равен нулю, а правее — от- 
рицателен и постоянен до ОПЗ. В ОПЗ поток ионов гидроксила так же, как поток ионов водорода, линейно рас- 
тет. В квазиравновесной ОПЗ поток ионов гидроксила резко меняется и удовлетворяет граничному условию 
(1) = „(т). Область рекомбинации является внутренним погранслоем для потоков водорода и гидроксила — 
здесь они резко меняются. 

Обсуждение и заключения. В данной работе рассмотрен перенос ионов бинарной соли через диффу- 
зионный слой у катионообменной мембраны. Предложена математическая модель названного процесса, кото- 
рая учитывает температурные эффекты, обусловленные реакциями диссоциации (рекомбинации) молекул воды 
и Джоулевым нагревом в растворе. Эта модель состоит из 9 обыкновенных дифференциальных уравнений с 
соответствующими краевыми условиями. Разработан алгоритм численного решения краевой задачи, приведено 


ее численное исследование. Установлены основные закономерности переноса ионов соли с учетом реакции 
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диссоциации (рекомбинации) молекул воды и температурных эффектов. 

Показано, что ВАХ неограниченно растет с увеличением скачка потенциала в запредельном режиме, 
причем значительно медленнее, чем в экспериментах. Следовательно, наряду с реакцией диссоциации (реком- 
бинации) воды существуют и другие механизмы сверхпредельного переноса (например, электроконвекция и 
гравитационная конвекция). Поэтому возникает проблема оценки взаимодействия реакций диссоциации (ре- 
комбинации), электроконвекции и гравитационной конвекции. 

Температурные эффекты от реакции диссоциации (рекомбинации) и Джоулева разогрева в ОЭН прак- 
тически незаметны (исключение — ОР). Джоулев нагрев в ОПЗ на два порядка больше охлаждающего эффекта 
реакции диссоциации воды. При рекомбинации в ОР выделяется примерно столько же тепла, сколько при Джо- 
улевом нагреве в расширенной ОПЗ. Однако из-за малых размеров ОР влияние этого тепла незаметно. Значит, 
можно считать, что есть только один источник тепла на межфазной границе в ОПЗ, который благодаря своему 
заметному размеру обусловливает значительное повышение температуры во всем диффузионном слое. Отсюда 
следует, что возможно возникновение и развитие гравитационной конвекции. 

В работе рассматривается раствор хлористого натрия, однако проведенное исследование справедливо 
для раствора любой бинарной соли. 

Общие выводы, следующие из полученных результатов, открывает возможность интенсификации про- 


цесса переноса ионов соли в электродиализных аппаратах. 
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